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Abstract — The paper deals with quantum pulse position mod- 
ulation (PPM), both in the absence (pure states) and in the 
presence (mixed states) of thermal noise, using the Glauber 
representation of coherent laser radiation. The objective is to 
find optimal (or suboptimal) measurement operators and to 
evaluate the corresponding error probability. For PPM, the 
correct formulation of quantum states is given by the tensorial 
product of m identical Hilbert spaces, where m is the PPM 
order. The presence of mixed states, due to thermal noise, 
generates an optimization problem involving matrices of huge 
dimensions, which already for 4-PPM, are of the order of 
ten thousand. To overcome this computational complexity, the 
currently available methods of quantum detection, which are 
based on explicit results, convex linear programming and square 
root measurement, are compared to find the computationally less 
expensive one. In this paper a fundamental role is played by the 
geometrically uniform symmetry of the quantum PPM format. 
The evaluation of error probability confirms the vast superiority 
of the quantum detection over its classical counterpart. 

Index Terms — Quantum detection, linear programming, 
square root measurement (SRM), least square measurement 
(LSM), geometrically uniform symmetry (GUS), thermal noise, 
pulse position modulation (PPM). 

I. Introduction 

Coherent optical communications represent the best way 
to convey an enormous amount of information over large 
distances, both through fibers and in free space. Quantum 
optical communications offer the possibility to improve the 
performance of (classical) optical communications by basing 
the design of the detection on the laws of Quantum Mechanics. 
Several authors have shown, both theoretically [3][16][24] 
and experimentally [4][5][19][23], that improvements, of at 
least a decade, can be obtained in terms of error probability. 
The uncertainty in quantum detection is intrinsically linked to 
the laws of quantum mechanics. In the language of classical 
optical systems, this uncertainty corresponds to the so called 
shot noise. Another source of uncertainty is the presence of 
thermal noise. 

Necessary and sufficient conditions for the optimal mea- 
surement set were found in pioneering papers by Holevo [6] 
and Yuen et al. [7]. Although the optimal measurement set 
is completely characterized, explicit solutions for the optimal 
measurement set are in general not available. It is therefore 
necessary to resort to a numerical evaluation based on convex 
semidefinite programming [20]. However, under a certain 
symmetry constraint, called geometrically uniform symmetry 
(GUS), a simple measurement, introduced by Hausladen et al. 
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[12], and known as square root measurement (SRM), turns 
out to be optimum. This SRM has the remarkable advantage 
that, starting from the given states, it can be straightforwardly 
evaluated. Moreover, even when it is not optimal, the SRM 
often provides "pretty good" upper bounds on the performance 
of optimal detectors. 

The problem of quantum detection in a noisy environment 
has received scarce attention in the literature, but the extension 
of SRM to mixed states [21] has opened new perspectives. 
In a previous paper [24], we applied this extension to the 
analysis of quantum communication systems, in the presence 
of thermal noise, using the Glauber theory on coherent states 
[1] to represent the noisy quantum channel. 

In this paper, the theory of [24], which was applied to 
PSK (phase shift keying) and QAM (quadrature amplitude 
modulation) systems, will be extended to PPM (quantum pulse 
position modulation) systems. The novelty of this modula- 
tion format is that it must be formulated (see [7] [19]) in 
a composite Hilbert space 3~C, given by the tensor product 
"Kf m of m identical Hilbert spaces Jfo, where m is the 
PPM's alphabet size. Thus, the quantum states correspond- 
ing to the classical binary PPM symbols are given by the 
tensor product of m Glauber coherent states. For instance, 
in 4-PPM, the composite states representing the symbols 
[0001], [0010], [0100], [1000], are given by 



|To) = |0)®|0)(gi|0)®|a) |7i) = |0}<8>|0)<8>|a}|<8>0) 
|T2> = |0}<8>|a}<8|0}<8|0) [73) = |a}<8> |0}<8> |0}<8> |0) 



(1) 



where the symbol 1 is represented by a Glauber state \a), 
with an average number of photons given by N s = |a| 2 , and 
the symbol by the Glauber "ground state" |0), which has 
zero photons. In the presence of noise (mixed states) the m 
composite states must be replaced by m composite density 
operators. If n is the dimension of the basic Hilbert space !Kq, 
the dimension of the composite space "K = 3if" Tl becomes 
N = n m , which gives rise to a problem of computational 
complexity. In the approximation of the Glauber representa- 
tion, in order to guarantee an adequate approximation, the 
choice of dimension n (theoretically infinite) depends on the 
average number of photons N s , which in turn depends on the 
range of the error probability P e that we wish to explore. To 
fix ideas, we anticipate the following schedule for 4-PPM 

• P e ~ 10~ 2 -> N s ~ 3 -> n ~ 10 , N = n m ~ 10 4 

• P e ~ 10~ 3 -> N, ~ 4.5 -> n ~ 15 , N = n m ~ 5 10 4 



10- 



N. ~ 6.5 



20 , N = n r 



16 10 4 



Thus, we realize that, while the dimension n of basic space 
can be confined to moderately small values, the dimension TV 
of composite space quickly becomes huge with consequently 
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severe numerical problems, mainly related to the eigendecom- 
position (EID) of TV x N matrices. Therefore, one of the main 
targets of the paper will be a systematic comparison of the 
currently available methods, to find the lowest computational 
complexity and, where possible, to search for analytical results 
for which the numerical evaluation does not represent a 
problem for any arbitrarily large dimension (see the EID of 
the PPM symmetry operator in Section III). 

The above schedule refers to a moderate amount of thermal 
noise (N = 0.05 thermal noise photons). Note that in the 
absence of noise, that is with pure states, the analysis is 
straightforward and the error probability is known, in closed 
form, from the pioneering works of Yuen, Kennedy and Lax 
[7], although the authors do not mention applications to PPM. 

This paper is organized as follows. In Section II we give 
a detailed overview of the quantum detection problem. In 
Section III we formulate the theory of PPM and prove that 
the quantum PPM format satisfies the GUS property, which 
surprisingly enough, has not been noticed elsewhere; the main 
task is to find an analytical expression for the symmetry 
operator and its EID. In Section IV we formulate the PPM 
quantum states according to the Glauber theory. Finally, in 
Section V, after a comparative discussion of computational 
complexity, we evaluate the PPM performances, up to the 
order m — 4. In Section VI we mention to the problem of 
implementation. 

II. Overview of Quantum Detection 

This overview is finalized to PPM. For a general overview, 
see [24]. 



A. Formulation of the problem 

The quantum detection problem is stated as follows [3]. A 
transmitter conveys classical information to a receiver through 
a quantum-mechanical channel, by choosing a quantum state 
from a set {pt\i = 0, . . . , m — 1} of density operators on an 
iV-dimensional Hilbert space "K. The density operators pi are 
chosen with given prior probabilities qt, but in what follows, 
we assume equal probabilities, i.e., qi = 1/m. The receiver 
detects the transmitted information with a set of positive oper- 
ator valued measurements (POVM) H Q , . . . , II TO _i. Provided 
that the state density operator is pi, the probability that the 
detection system reveals the state j, is given by the transition 
probability p(j\i) = Tr^IL,), i,j = Q,...,m—l. Then, the 
probability of correct detection becomes 

1 rn — 1 1 m — 1 

i=0 i=0 

The description through density operators represents the gen- 
eral case, and includes the pure state case, in which pt reduces 
to the rank-one operator p, = |7»)(7i|. 

Quantum detection is simplified and has specific properties 
if the state constellation exhibits the Geometrically Uniform 
Symmetry (GUS), that is if a unitary operator S exists with 



the properties 1 

Pl = S*p S-* , S m = I K (3) 

where S is the symmetry operator, and po is the generating 
density operator. The simplifications with GUS are: 1) The 
quantum source is completely specified by po and S. 2) 
The POVM Hi can be chosen to have the same GUS as 
the density operators pi [21], namely, IL; = S l U S^\ The 
search can then be limited to the reference POVM IIo. 3) The 
transition probabilities p(j\i) depend only on the difference 
i — j, mod to (their matrix is circulant). In fact, p(j\i) = 
Tr [n St-i po S-^-^] . Consequently 

P c = Tr( Po n ) . (4) 

The positive semidefinite (PSD) Hermitian operators p 
and il may, in some respect, be redundant and can then 
be investigated through their factors. In an iV-dimensional 
Hilbert space the density operators po are represented by 
N x N complex matrices and can be factorized in the form 
Po = 7o7o> where 70 is a convenient N x ro matrix, with 
ro = rank(po)- This factorization is not unique, but the 
ambiguity turns out to be irrelevant in quantum detection [21]. 
In particular, if the states are pure, then r = 1 and 70 reduces 
to the ket I70). Since the rank of the optimal POVM IIo is 
not greater than ro [21], we postulate a similar factorization 
IIo = Mo Mo °f me measurement operator. In the presence 
of GUS the factors of pi and IL can be obtained from the 
reference factors as 7$ = S l 70 and pn = S l po, where 7$ and 
Pi have the same dimension N x r . 

B. Optimal solution (minimum error probability) 

In their pioneering papers, Holevo [6] and Yuen et al. [7] 
found necessary and sufficient conditions (see also [21]) for an 
optimum POVM set (which minimized the error probability). 
In a few special cases, the optimal POVM set is explicitly 
known, namely: 1) for binary systems in the general case of 
mixed states, and 2) for m-ary systems with pure states having 
GUS. Case 2) will be seen below in the context of SRM. Case 
1) is due to Helstrom [3], who obtained the following explicit 
result: 

Proposition 1 Let p and p\ be the density operators of 
a binary quantum system. Define the decision operator D = 
QiPi — qoPo, where qo and qi are the prior probabilities, and find 
the corresponding EID: D = |ei) (e^ |, where the eigen- 

values 6i are real and |e») are the corresponding eigenvectors. 
Then, the optimal POVM are given by YLq = J2e <o \ £ i)( e i\ 
and ili = ]C e >o l e «)( £ il- The corresponding correct detection 
probability is given by P c = qo + J2 u >o e i- 

In the general case, explicit solutions for the optimal POVM 
are not available. However, the problem of maximizing the 
correct detection probability P c can be solved numerically by 
convex semidefinite programming [20]. For instance, one can 
use the LMI (Linear Matrix Inequality) Toolbox of MatLab, 

'in this formulation equal prior probabilities are implied. In the general 
case, the weighted density operators qiPi should be considered [21]. Note 
that if the GUS holds with equal probabilities, as in (3), it does not hold with 
arbitrary probabilities. 
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which permits evaluating the optimal performance with any 
desired accuracy. Of course, if the dimensions are very large, 
as in the case of PPM with mixed states, the numerical 
evaluation may become very time and storage consuming. 

C. SRM solution (square root measurement) 

A more straightforward, although not optimal, approach to 
the problem is offered by the SRM. The approach for pure 
states was proposed by Hausladen et al. [12] and thoroughly 
investigated by Eldar and Forney [18]. The generalization to 
mixed states is due to Eldar et al. [21]. In the last paper 
it is shown that SRM is equivalent to LSM (least square 
measurement). The SRM formulation is stated in terms of the 
state matrix T — [70, 71, ... , 7 m -i ] and of the measurement 
matrix M — [/j,q, p,\, . . . , p, m -i ], both of dimension Nxmra. 
Now, the quantum detection problem can be reformulated in 
terms of the matrix T (given) and the matrix M (to be found). 

Proposition 2 (see [18]) In the SRM the measurement 
matrix M is given by the two equivalent expressions 



M = r 1/2 r, m = tg- 1/2 



(5) 



where T x l 2 and G 1 / 2 are the inverse square roots (in the 
Moore-Penrose generalized sense) ofT and G, respectively. 

The alternative approaches to evaluate M, either via T -1 / 2 
or via G -1 / 2 , are important for efficient computation. We 
develop these results in the case of GUS, which is of special 
interest for PPM. The first of (5) gives the reference POVM 
as n = T- 1 ' 2 p^T- 1 / 2 . Once T" 1 / 2 is evaluated, we 
can obtain the transition probabilities, and then the correct 
detection probability from (4), which becomes 



P c = Tr UpoT- 1 / 2 ) 2 



(6) 



The Gram matrix approach is less direct. The matrix G is 
block circulant and its decomposition is related to the discrete 
Fourier transform (DFT). In [24] we found: 

Proposition 3 In the presence of GUS thei,j block of the 
Gram matrix, of size r x r , can be written in the form 



SRM optimality with GUS. Under certain conditions, SRM 
provides the optimal solution for the general case of mixed 
states. A (sufficient) condition is that the state constellation 
has GUS and the reference factors satisfy the relation [21] 
Po 7o = a Ih > where a is an arbitrary constant. Unfortunately, 
for PPM this condition holds only for pure states. 

III. Quantum Pulse Position Modulation (QPPM) 

QPPM is an attractive modulation format, which has recently 
been considered [22] for possible applications to deep space 
communications. Indeed, it appears well suited to existing 
laser modulation techniques, in that it can be regarded as a 
coded version of the OOK format. In particular the require- 
ments for the transmitting laser are the same as for the OOK 
coherent modulation, with lower average power. On the other 
hand, the model of the QPPM we consider below is a paradigm 
for a whole class of quantum communication scheme (as 
wavelength modulation). Moreover, it has been proved [3] 
that the model is a quantum-mechanical equivalent of a 
general communication system with (classically) orthogonal 
waveforms. 

The performances of QPPM systems are known only in 
very particular cases. The minimum error probability for 
QPPM with quantum detection in the absence of thermal noise 
detection (see below) was found in 1975 by Yuen et al.. In 
1976 Helstrom [8, Chapter 6] evaluated the error probability 
in the presence of thermal noise, by assuming that a classical 
OOK receiver counts the photons slot-by-slot and decides 
according to a majority rule. In 1983 Dolinar [5] suggested an 
adaptive measurement approach, called conditionally nulling, 
which achieves very nearly the same error probability of the 
optimum quantum receiver in the absence of thermal noise. 
However, no precise evaluation exists in the literature for 
optimal performances of PPM in the presence of thermal 
noise. 

In QPPM the quantum-mechanical channel is modeled by a 
composite Hilbert space "K, given by the tensor product Mf 1 ™ 
of to identical Hilbert spaces "K of dimension n. In this space 
the transmitter chooses one of to pure quantum states, given 
by the tensor products 



Gij = 7*7, = 7o^'"So = ^ E WW-»E k (7) hi) = hi,m-i)( 



l7ii)®l7io) 



0. 



■1 (10) 



fc=0 



where W m — e j27 7 m and the matrices E k , of order H, have 
the alternative expressions 

m—l 

E k = to 7o * Y k 70 - ]T G M W~ k * . (8) 



i=0 



Here, Y k are the matrices appearing in the EID of S (see (14)). 
The correct detection probability is finally given by 



where = I7 ) for j ^ i and ^y) = I7 1 ) for j = i. 
In a practical QPPM scheme the two states are assumed to 
be I7 ) = |0) (the Glauber ground state) and I7 1 ) = |a) (a 
Glauber state with an average number of photons N s = |a| 2 ). 
The explicit constellation for m — 4 was given in (1). If the 
states are mixed, the composite pure states must be replaced 
by the composite density operators 



Pi = Pi,m-1 ® • • • ® Pn ® PiO 



i = 0, 



,m-l (11) 



1 m— 1 „ 



(9) 



fe=0 



In conclusion, with the Gram matrix approach, the perfor- 



mance evaluation requires finding the square roots E k 
the EID, of the ro x ro matrices E k defined by (8). 



1/2 



, by 



where pij = p° for j ^ i and pij = p 1 for j = i, and p° and 
p 1 are assigned density operators on the space J£o. 

A. Symmetry Operator 

The GUS structure of the QPPM constellation is apparent. 
Naively, the state p i+i is obtained from p, through a symmetry 
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operator S that pushes each Kronecker factor by one position 
(modulo m) to the left. However, to translate the above words 
into a formula is not a trivial problem, since the operator S 
is not separable into a tensor product of operators. 

Proposition 4 The density operators pi of QPPM form 
a GUS constellation, namely pi — S l poS~ l . The symmetry 
operator is given by 

n-1 

S=J2 w n(k)®I H ®<(k), (12) 

fc=0 

where w n (k) is the column vector of length n, whose elements 
are with exception of a 1 in the k-th position, and Ih is the 
identity matrix of order H = n m_1 . 

As an illustrative example, the symmetry operator for n = 2 
and m — 4 (4-PPM) takes the form 



rl 
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PROOF. The result is a particular case of a property of the 
Kronecker product (see [9]). Given two matrices A and B, 
with dimensions r x h and s x k, respectively, the Kronecker 
products A (g) B and B (g) A are different, although they have 
the same dimension rs x hk and contain the same entries, in 
different positions. Moreover, the products are related by 

B® A= S rs [A®B]S-l , (13) 

where the matrices S rs and Shk are permutation matrices 
of order rs and hk, respectively, known as perfect shuffle 
matrices, which are independent of the particular matrices A 
and B and depend only on their dimensions. The matrix S ra is 
given by S rs = J2k=o w s (k)(g>I r (g>w'*(k). Now, the symmetry 
operator of QPPM gives the mapping 

io] * [ Pi,m— 20 * • -®piQ ] ®Pi,m— 1 

which is a particular case of the general transformation (13) 
with square matrices A = Pi, m -\ of order n and B = 
Pi.m-2® ■ ■ -®Pio °f order H = n m_1 . By putting r = h = n 
an s = k = H into (13), we have S rs = Shk = S n H and we 
find the expression (12). 

B. Eigendecomposition of the Symmetry Operator 

As shown above, in the performance evaluation of a quan- 
tum detection system that satisfies a GUS constraint, the EID 
of the symmetry operator S is a key step. Since S is unitary 
and S m = /, the EID has the form 

m— 1 

S=Y J ^Y*Yk (14) 

fc=0 

where the Y k collect, as their columns, the kets \yf) corre- 
sponding to the eigenvalue = W~ k . The problem is the 



evaluation of the matrices Yk- We can also exploit the fact that 
in our case S reduces to a permutation matrix. The EID of 
permutation matrices is an intricate topic (see for instance [2]), 
at least for the general case. However, the particular structure 
(12) of the operator simplifies the results. In any case, the 
cycles generated by S, according to the following definition, 
play a fundamental role for the EID. Provided that w(k), 
k = 0, . . . , n m — 1 are the unit vectors of the space % — %® n , 
the distinct unit vectors w(ko),w(k\), . . . ,w(k p -\) form a 
cycle (of period p) if S p w(ko) — w(fc ) and, for each 
i = 0, . . . ,p — 1, w(ki) — S l w(k ) . 

In Appendix A we discuss the properties of the symmetry 
operator S listed in the following proposition. 

Proposition 5 1) The periods of the cycles of S are divisors 
of m. 2) The indexes k n , . . . , fc r „_i that identify the cycle 
w(k n ), w(ki), . . . , w(k p -i) are obtained from k n through the 
recursive relation 

k i+1 = nh mod (n m - 1) . (15) 
3) For any cycle w(k ) , w(ki) , . . . , w(fc p _i) the vector 

1 

Vj = —J2 W i hw (kh) j = 0,...,p-l (16) 
vP h=0 

is an eigenvector of S associated to the eigenvalue W~K The 
eigenvectors so obtained constitute an orthonormal basis of the 
space "K. 

As a conclusion, we have the following EID of the symme- 
try matrix S 

S=EE Y^e-^tvV c] V; 3 (17) 

where the subscript p \ m denotes that the sum is extended to 
the divisors of m (m included), % p is the set of cycles of 
minimum period p and, for each cycle c e X p , the vectors 
V c j, j = 0, . . . ,p — 1 are the eigenvectors generated by the 
cycle c. After suitably reordering indexes, (17) can be written 
in the form 

m— 1 

s=Yl W ^ Y * ( 18 ) 

i=0 

where YjY* turns out to be the projector mapping the space 
IK into the autospace associated to the eigenvalue W^. 

In order to evaluate the EID (17), we can use the following 
proposition, proved in Appendix B. 

Proposition 6 The number N p of unit vectors of minimum 
period p is given by the recursive relation n p — J2d\ p Nd- 
The number of cycles of period p is N p /p. The multiplicity 
of the eigenvalue W~ h is given by n h = E(m/ P ) \h N p/P' 
h = 1, 2, . . . , m, where the summation is extended to all the 
p such thatm/p divides h, and n m gives n . 

C. Examples 

We reconsider the previous example of S, for n = 2 and 
m = 4 (4-PPM). The possible periods for the cycles are 1, 2, 
and 4. By applying the recursive relation (15), which becomes 
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= 2ki mod 15, we find the cycles {0} and {15} of 
period 1, {5,10} of period 2, {1,2,4,8}, {3,6,12,9}, and 
{7, 14,13, 11} of period 4. 

The eigenvalues 1, W 4 , W%, Wf have 6, 3, 4, and 3 eigen- 
vectors, respectively. The matrices collecting the eigenvectors 
associated to the eigenvalues 1, W4, Wf , Wf are then 



of a single radiation mode is given by the ket 



v5 




Y 1 = 



Y 2 



V2 






1 

^2 



Y 3 = 



As another example, let n = 2 and m = 10, so that n m — 
1024. The possible periods for the cycles are p = 1, 2, 5, 10. 
The number of unit vectors with period p are Ni = 2, N 2 = 2, 
N 5 = 30, and N w = 990. The numbers N p /p are 2, 1, 6, 99, 
respectively. The matrices have dimension 1024 xn„ where 

n = 108, Til = "3 = "7 = n 9 — 99, «2 = n A — n 6 = n 8 = 

105, and n 5 = 100. 

It can be shown that the sizes of the matrices involved 
increase very rapidly. Fortunately, the matrices Yi are very 
sparse and strongly structured. Several appropriate tricks can 
then be exploited to alleviate the computational complexity. 

IV. Quantum states in coherent optical 

COMMUNICATIONS 

In quantum optical communication systems the correct setting 
of quantum states is provided by the Glauber's celebrated 
theory [1]. In the case of PPM, this theory will provide the 
basic density operators p° and p 1 , representing the symbols 
and 1, from which the m density operators pi can be obtained 
by tensorial products. 



A. Density Operators from Glauber's Theory 

In this theory, a coherent state, representing a monochro- 
matic electromagnetic radiation produced by a laser, is for- 
mulated in an infinite dimensional Hilbert space "K equipped 
with an orthonormal basis \n),n — 0, 1, 2, . . ., where \n) are 
called number eigenstates. Each state \n) is said to contain 
exactly n photons. In this context the Glauber representation 



|a> = e-*H 2 £^|n) 
n vn! 

n— 



(19) 



where a is the complex envelope that specifies the mode and 
N a — \a\ 2 represents the average number of photons, when 
the system is in the coherent state |a). 

The representation (19) is valid when the receiver observes 
a pure state with a known parameter a, which in the context 
of communications may be regarded as the signal. In the 
presence of thermal (or background) noise, the signal becomes 
uncertain and must be represented through a density operator. 
In this case the Glauber theory provides an infinite matrix 
representation ||p mn || of the density operator, given (for < 
m < n) by 



M = (i 




-(l-v)\ a \ 2 



(20) 



K(K+1) / 

where N represents the average number of photons associated 
with the thermal noise, v = N/(l + N), and L™~ m (x) are 
generalized Laguerre polynomials. The entries for m > n 
are obtained by p nm {a) = p* mn {a). The diagonal elements 
p mm (a) obey a Laguerre distribution and give the probabil- 
ities of exactly m photons being present, when the quantum 
system is in the noisy state represented by the density operator 
p(a) [3]. For the ground state I7) = |0) the above expression 
degenerates and the corresponding matrix representation be- 
comes diagonal, namely p mn (0) = S mn (l — v)v n . 

The infinite dimensional matrix | |/o mn (a)| |, < m,n < 
00 gives a correct representation of the density operator, but, 
for practical calculations, we introduce a finite dimensional 
approximation by truncating to n £ terms. For the choice of 
n e , to achieve a given accuracy, we follow the quasi-unitary 
trace criterion, proposed in [24], which is based on the fact 
that a density operator has unitary trace. Then, we choose n € 
as the smallest integer, such that El=o Pramftt) > 1 — e, 
where e is the required accuracy. Thus, for a given e, n e can 
be evaluated using the Laguerre distribution p mm (a). 

Numerical accuracy. The above numerical approach implies 
a twofold approximation. First, the truncation of the infinite 
matrix to a finite n £ x n £ matrix p. It is reasonable to 
think that, at the increasing of n £ , the approximate results 
converges to the correct ones. However, the evaluation of the 
corresponding error seems to be a difficult problem and we 
made no attempt to evaluate it. A second approximation is 
that the truncated density matrix is not normalized to unitary 
trace. Of course, it can be normalized as p = cp, where 
c = 1/Tr(p) < 1/(1 - e). In the PPM the truncation is 
operated on the basic operators p° and p 1 , which we normalize 
as p° = cp° and p 1 = cp 1 . To see the non normalization 
effect on the error probability P e = 1 — P c , we consider the 
SRM approach stated by Proposition3, and mark with a bar the 
quantities obtained from normalized density operators. For the 
reference factor we get 70 = c m / 2 70 and for the square-root 
matrices E l k /2 = c m / 2 E^ 2 , and finally from (9) P c = c m P c . 
Hence, P e = (1 — c m ) + c m P c and, considering that e«l 
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and P e <C 1, we find P e — P e ~ me. In conclusion, if we 
choose e <C P e , we see that normalization is irrelevant. In the 
numerical evaluation of Section V, where P e > 10~ 6 , we have 
chosen e = 1CT 8 . 

B. Factorization of the density operators. Choice of the rank 

Once the finite n e x n € approximation of the density operator 
has been established, we need a factorization of the form 
p(a) — 7(01)7* (a) for a suitable matrix 7(a). For a = 
(ground state) the factorization is immediate, since p(0) is 
diagonal, and we find 7(0) = \J p(0) = \\S mn ^(l — v)v n \\. 
For a / we evaluate the reduced EID of p(a), namely 
p{a) = ZfiAlZ^, where h is the rank of p(a), Zh is an 
n e x h matrix that collects the eigenvectors corresponding to 
the h positive eigenvalues A 2 , and is the h x h diagonal 
matrix that collects the Af . Hence, 7(a) = Z^Ah is a factor 
of p(a). 

A critical step in the numerical evaluation is the choice 
of the rank h, given by the number of numerically relevant 
eigenvalues. To clarify the problem we develop a specific case: 

a = V5 , N s = 5 , N = 0.1 , e = 10~ 5 -» n e = 20. 
Now, in theory, p(a) has full rank h — n e , as can be seen from 
the list of its eigenvalues, obtained with a high accuracy 

0.150285 , 0.00231095 , 0.0000353779 , 5.20725 10~ 7 
7.24874 10 -9 , 9.45157 10 -11 , 1.14603 10 -12 
3.10867 10~ 21 , 2.32186 10~ 21 , 1.37631 10" 



3.01775 10~ 25 



but in practice, we can limit the list to the first 3 eigenvalues, 
neglecting the rest, which means that we assume h = 3 as 
a "practical" rank. As a check, the reconstruction of p(a) 
obtained in this way assures an accuracy < 10~ 8 . To find 
the "practical" rank in the general case, we consider the 
reconstruction error Ap = p — jhJh> wnere the factor is 
obtained by taking only h eigenvalues. Then, we can evaluate 
the maximum error, or the mean square error, as a function 
of h and choose h = h v to achieve a given accuracy v (for a 
detailed discussion of n e and h v versus the accuracies e and 
v, see [24]). 

V. Application examples 

We apply the theory of the previous sections to evaluate the 
error probability in QPPM systems. But first we compare the 
developed methods to give a guide to possible ranges in the 
different numerical evaluations. 

A. Comparison of computational methods 

In the previous section, we saw that, given the accuracies 
e and v, we can evaluate the basic density operators p° = 
P(0) = 7h(0)7h(0)* and p 1 = p{a) = Jh{a)lh{a)* for any 
choice of signal parameter N s = \a\ 2 and of thermal noise 
parameter N (it is expedient to choose the same "practical" 
rank h both for p° and p 1 ). Hence, the dimensions and the 
ranks are as follows: 

• 7 and 7 1 n x h rank h, p° and p 1 nxn rank h, 

• 7i N x H rank H, p { N x N rank H, with 
N ~n m and H — h m , 

• state matrix T N x mH rank r — mH (if mH < N), 

• Gram operator T N x N rank mH, Gram matrix G 
mH x mH rank mH. 



Having this schedule in mind, we now compare the different 
evaluation methods, considering that the practical limit is 
determined by the dimensions of the input data in convex 
linear programming (CLP), and that in general, it is given by 
the dimensions of the matrices to be eigen-decomposed. 

1 ) Optimal solution: This solution is available for: 

• pure states (N = 0), for any order m, 

• mixed states, for to = 2 from Helstrom's theory, but an 
N x N EID is required, 

• CLP with input the N x N density operators p,. 

2) SRM solution: This solution is always available, for any 
order to, for both pure and mixed states. We have two possible 
approaches 

• via Gram operator, requires the EID of an N x N matrix, 

• via Gram matrix, requires the EIDs of H x H matrices. 

Since H <C N, we have used the Gram matrix approach, with 
a practical limit for the EID computation of H ~ 1500. With 
the CLP (Matlab toolbox LMI), we found the practical limit 
of the size N of the input density operators to be no greater 
than 150. These limits refer to a standard personal computer. 
To complete the preliminaries we consider two more topics. 

Error probability with pure states. We saw that SRM gives 
the optimal solution with pure states, by virtue of the PPM's 
GUS. For the explicit evaluation, we use (8), where Gq s = 

7o7s = (7ol7s> = (7oo|7«o) ■ ■■ (7o,m-i|7s,m-i)- Recalling 
that the inner product of two Glauber states is given by [1] as 
(a\/3) = exp(-i(H 2 + |/3| 2 - 2a* /3)), we get G 0s = 1 for 
s = and Go s = exp(— \a\ 2 ) = exp(—N s ) for s^0. Hence, 
from (9), we have 

2 



TO 2 (V 



1 + (to - 1) e"^ + (to - l)Vl~- 



-n 3 



(21) 

which is in perfect agreement with the result of [7], obtained 
with a different approach. 

Comparison with classical optical PPM. In a classical 
optical PPM system based on a photon counter, the error 
probability was evaluated by Helstrom [8, Chapter VI], using 
the Laguerre distribution of photon arrivals. He found 



lit 



exp 



(1- V )(1 



l-v l 



where v = K/(l + K). 



B. Application to 2-PPM 

For to = 2, we do not encounter numerical problems 
because, for the range of interest, the dimension are relatively 
small. We first present the results in the absence of thermal 
noise (pure states), compared with the performance of a classic 
optical PPM system (Fig. 1), which confirms the substantial 
improvement of the quantum system with respect to the 
classical one. 

In the presence of thermal noise, for 2-PPM (as for any 
other binary format) an exact evaluation is possible using 
Helstrom's theory (see Proposition 1). This possibility was 
also used to check the results obtained with the Matlab LMI 
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Fig. 1 



Error probability of classic and quantum 2-PPM in the 
absence of thermal noise. N a is the average number 
of photons per symbol. 
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Fig. 2 



Error probability in 2-PPM versus the average number 
of photons per symbol N a for some values of the ther- 
mal noise parameter N. Comparison of SRM results 
with Helstrom's optimum results. 
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Fig. 3 



Error probability versus the average number of pho- 
tons per symbol N„ of classic and quantum 3-PPM 
in the absence of thermal noise. 
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Fig. 4 — Error probability in 3-PPM versus the average number 
of photons per symbol N s for some values of the ther- 
mal noise parameter X. Solid lines refer to quantum 
detection and dashed lines to classical photon counter 
detection. 



toolbox. For SRM we used the Gram matrix approach. The 
results are shown in Fig. 2 for some values of the thermal 
noise parameter N and show that SRM overestimates the error 
probability by about 30%, with respect to the optimum. We 
conclude that SRM is also "pretty good" in the presence of 
thermal noise. 

C. Application to 3-PPM 

For m = 3, we encounter some numerical problems that 
can be overcome using the alternative approaches. We first 
present our results in the absence of thermal noise (pure 
states), compared with the performance of a classic optical 
PPM system (Fig. 3). We realize the considerable improvement 
of the quantum system with respect to the classical one. 

In the presence of noise, we used the SRM Gram matrix 
approach. The results are shown in Fig. 4 as a function of N s , 



for some values of the thermal noise parameter N. To extend 
the error probability range down to 10~ 5 we chose n = 40, 
h = 8, hence H = 8 3 = 512, which is within the limit of our 
computational capability. The figure also compares quantum 
detection with classical photon counter detection, and shows 
more than 3 decades improvement in error probability. 

Since SRM is not optimal in the presence of noise, for the 
3-PPM we tried to evaluate the minimum error probability by 
CLP. As shown in Fig. 5, this was feasible for a very small 
range of N s , where the dimension N ~ n 3 can be kept below 
a value of 150. Anyway, this has permitted to show that the 
SRM gives a very small overestimate of the error probability. 

D. Application to 4-PPM 

For m = 4, the numerical problems become very severe 
and the comparison with the optimal case (via CLP) is not 
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Fig. 5 — Error probability in 3-PPM versus the average number 
of photons per symbol N s for the value X = 0.1 
of the thermal noise parameter. Solid line refers to 
SRM computation and dotted lines to convex linear 
programming (CLP) computation. 
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detection and dashed lines to classical photon counter 
detection. 
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possible. We present the results in the absence of thermal 
noise (pure states), compared with the performance of a classic 
optical PPM system in Fig. 6 and we realize the considerable 
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Fig. 6 — Error probability versus the average number of pho- 
tons per symbol N s of classic and quantum 4-PPM 
in the absence of thermal noise. 



improvement of quantum 4-PPM with respect to the classical 
one. 

In the presence of noise, we used the SRM Gram matrix 
approach. The results are shown in Fig. 7 as a function of N a 
for some values of the thermal noise parameter N. To cover 
the error probability range down to 10~ 5 we chose n = 30, 
h = 6, hence H = 6 4 = 1296, which is at the limit of 
our computation capability. The figure also compares quantum 
detection with classical photon counter detection and shows 
more than 3 decades improvement in error probability. 



VI. On the Implementation 

The problem of designing and implementing an optimal 
quantum receiver is a very difficult one also in the simplest 
cases, as witnessed by the history of the optimal detector 
for the quantum OOK modulation format. As known, the 
error probability that can be obtained with a classical detector 
(direct photon counter) in the absence of thermal noise is 
given P e = (l/2)e _2JVs , while the error probability achiev- 
able by optimal quantum detection (the Helstrom bound) is 
(1/2) [l - x/T^F 
(l/4)e- 2Ar = 



-2N S 



with asymptotical approxi- 



opt 

mation P op 

In 1973 Dolinar [5] proposed an optimum receiver ideally 
achieving the Helstrom bound. Under the Dolinar's receiver 
it was the idea of using an adaptive detection, by adding to 
the received coherent field a time-varying feedback controlled 
field. Owing to the lack of an adequate technology in fast 
electro-optics, the realization of the Dolinar receiver was not 
possible, until in 2007 Cook et al. [23] demonstrated its 
feasibility in a significant range of physical parameters. 

For a quantum PPM a natural solution could be to use 
an OOK receiver for a slot-by-slot hard detection slot by 
slot. For pure quantum states |0) and \a), a classical photon 
counter receiver may fail only in detecting the state \a). It 
is easy to show that the error probability turns out to be 
P e = (m/(m — l))e~ Ns . This is an exponentially degraded 
performance with respect to the quantum limit (21) having 
asymptotical approximation P opl w (1/4) (in — l)e~ 2Ns . 

On the other hand, an algorithm devised by Dolinar and 
called "conditionally nulling" [10] allows to obtain an asymp- 
totical error probability P op , w (l/2)(m — l)e~ 2Ns , only 
double of the quantum limit. This near-optimum algorithm 
decides at each slot (on the basis of the previous results) 
whether or not to add to the incoming optical signal a negative 
replica of the signal corresponding to the quantum state \a). 

Unfortunately, the above considerations hold true only in 
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the absence of thermal noise. The degradation due to thermal 
noise is object of research in progress. 

VII. Conclusion 

We have examined the quantum PPM format, both in the 
absence (pure states) and in the presence (mixed states) of ther- 
mal noise, using the Glauber representation of coherent states. 
Our target was a numerical evaluation of the performance of 
4— PPM, where the dimension of the Hilbert space reaches ten 
thousand and more. To this end, we explored all the currently 
available methods, namely: exact solution, convex linear pro- 
gramming (CLP), and square root measurement (SRM), and 
also exploited the geometrically uniform symmetry (GUS), a 
property not considered before for the PPM, with the goal to 
find the least expensive implementable solution. We succeeded 
in this, but just at the limit, so that the methods used do not 
appear to be useful for higher orders, e.g., 8-PPM. The use of 
CLP was limited to a small error probability range in 3-PPM. 

Whenever possible, we compared SRM with optimal de- 
tection and concluded that SRM is very close to optimality. 
Comparisons with the performance of classical detection con- 
firms the superiority of the quantum detection also when the 
thermal noise is present. 

Appendix A: Proof of Proposition 5 

Given a unit vector w(k) of the composite Hilbert space 
'K = 5Cf"\ a positive integer q is a period of w(k) (with 
reference to the operator S) if S q w(k) = w(k). The minimum 
p of periods of w(k) is its minimum period. Simple consid- 
erations guarantee some obvious facts, such as the existence 
of a minimum period for each unit vector w(k) and the fact 
that S q w(k) = w(k), if and only if, p is a divisor of q. In 
particular, since S m w(k) = w(k), the minimum period p of 
each unit vector must be a divisor of to. 

If we use the n-ary representation of the integer k, namely, 
k = fc m _in m_1 + . . . + kin + k , with < h < m, the 
permutation 



S[w n {k m -{) <g> . . . <g> w n (ki) <g> w n (k )] 

= w„(fc m _ 2 ) ® • • • ® w n (k ) <g> w n (k m -i) 



(22) 



can be written in the form Sw(k) = w(h), where h = 
k m -in m ~ x + . . . + k n + fc TO _i. Since 

m-2 

nk = A: i n i+1 +fc m _ 1 +/c m _ 1 (n m -l) = h mod (n m -1) , 

the recursion (15) holds. 

Since the vectors (16) satisfy the condition 



SVi = 



1 



p-i 



:'^^^ h/P Sw{k h ) 
V P h-n 



h=0 
p-2 

£< 

h=0 



1 

e -i2nj/py 



,i2njh/p 



W(kh+l) + e 



i2nj(p-l)/p 



tu(fco) 



Vj is an eigenvector of S associated to the eigenvalue W p ■? = 
e~ l2v ^ p . Since w^.Wkj = Sij, it can be immediately verified 



that the eigenvectors corresponding to a single cycle are or- 
thonormal. Moreover, eigenvectors corresponding to different 
cycles are orthogonal, in that they are linear combinations of 
disjoint sets of unit vectors. Then, item 3) follows. 

Appendix B: Proof of Proposition 6 

By examining the permutation (22), it turns out that the unit 
vector w(k) has period d (not necessarily minimum), if and 
only if, the n-ary kok\ . . . fc TO _i representation of the index 
k is periodic, of period d, with respect to the cyclic shift, i.e., 
it is given by the concatenation of m/d equal subsequences 
fcofci . . . kd-i- Since the number of such subsequences is n d , 
this is also the number of unit vectors with period d. Since a 
unit vector has period d, if and only if, its minimum period 
divides d, the recursion n p = ^ d | p Nd follows. Moreover, 
for each subsequence k ki . . . kd-i, its d cyclic shift give rise 
to the same cycle, so that the number of cycles of period d 
is N d /d. Finally, the multiplicity is stated from the fact that 
a cycle of period p gives rise to an eigenvector associated 
to W^ h , if and only if, there exists an integer r such that 
yy-rm/p _ w~h^ j e 5 jf an( j on jy jf^ /j, i s a multiple of m/p. 
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